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Abstract 

Iterative rounding and relaxation have arguably become the method of choice in dealing with 
unconstrained and constrained network design problems. In this paper we extend the scope of the 
iterative relaxation method in two directions: (1) by handling more complex degree constraints in the 
minimum spanning tree problem (namely laminar crossing spanning tree), and (2) by incorporating 
'degree bounds' in other combinatorial optimization problems such as matroid intersection and lat- 
tice polyhedra. We give new or improved approximation algorithms, hardness results, and integrality 
gaps for these problems. 

• Our main result is a (1, 6+0(log n)) -approximation algorithm for the minimum crossing span- 
ning tree (MCST) problem with laminar degree constraints. The laminar MCST problem is a 
natural generalization of the well-studied bounded-degree MST, and is a special case of gen- 
eral crossing spanning tree. We also give an additive 17(log" m) hardness of approximation for 
general MCST, even in the absence of costs (a > is a fixed constant, and m is the number of 
degree constraints). 

• We then consider the crossing contra-polymatroid intersection problem and obtain a (2, 26 + 
A — 1) -approximation algorithm, where A is the maximum element frequency. This models for 
example the degree-bounded spanning-set intersection in two matroids. Finally, we introduce 
the crossing lattice polyhedra problem, and obtain a (1,6 + 2A — 1) approximation under 
certain condition. This result provides a unified framework and common generalization of 
various problems studied previously, such as degree bounded matroids. 

1 Introduction 



Iterative rounding and relaxation have arguably become the method of choice in dealing with uncon- 
strained and constrained network design problems. Starting with Jain's elegant iterative rounding scheme 
for the generalized Steiner network problem in [18], an extension of this technique (iterative relaxation) 
has more recently lead to breakthrough results in the area of constrained network design, where a number 
of linear constraints are added to a classical network design problem. Such constraints arise naturally in 
a wide variety of practical applications, and model limitations in processing power, bandwidth or budget. 
The design of powerful techniques to deal with these problems is therefore an important goal. 

The most widely studied constrained network design problem is the minimum-cost degree-bounded 
spanning tree problem. In an instance of this problem, we are given an undirected graph, non-negative 

*IBM T.J. Watson Research Center. {nikhil,rohitk:, viswanath} @us.ibm.com 
^University of Waterloo. jochen@uwaterloo.ca 
■''Technische Universitat Berlin, peis@math.tu-berlin.de 



1 



costs for the edges, and positive, integral degree-bounds for each of the nodes. The problem is easily 
seen to be NP-hard, even in the absence of edge-costs, since finding a spanning tree with maximum 
degree two is equivalent to finding a Hamiltonian Path. A variety of techniques have been applied to 
this problem IIIIISIIBIEIIEZIIIIIIIII, culminating in Singh and Lau's breakthrough result in ffH- They 
presented an algorithm that computes a spanning tree of at most optimum cost whose degree at each 
vertex v exceeds its bound by at most 1, using the iterative relaxation framework developed in 112411311 . 

The iterative relaxation technique has been appUed to several constrained network design problems: 
spanning tree lOTI . survivable network design |[24l l25l . directed graphs with intersecting and crossing 
super-modular connectivity Il24l l4l. It has also been applied to degree bounded versions of matroids and 
submodular flow |[T9l . 

In this paper we further extend the applicability of iterative relaxation, and obtain new or improved bi- 
criteria approximation results for minimum crossing spanning tree (MCST), crossing contra-polymatroid 
intersection, and crossing lattice polyhedra. We also provide some hardness results and integrality gaps 
for these problems. 

Notation. As is usual, when dealing with an undirected graph G = {V,E), for any S C 1/ we let 
6g{S) := {{u,v) ^ E \ u ^ S, V ^ S}. When the graph is clear from context, the subscript is dropped. 
A collection {Ui, • • • , Ut} of vertex-sets is called laminar if for every pair Ui, Uj in this collection, we 
have Ui C Uj, Uj C Ui, or Ui n Uj = 0. A {p, f{b)) approximation for minimum cost degree bounded 
problems refers to a solution that (1) has cost at most p times the optimum that satisfies the degree 
bounds, and (2) satisfies the relaxed degree constraints in which a bound b is replaced with a bound f{h). 

1.1 Our Results, Techniques and Paper Outline 

Laminar MCST. Our main result is for a natural generalization of bounded-degree MST (called Lam- 
inar Minimum Crossing Spanning Tree or laminar MCST), where we are given an edge-weighted undi- 
rected graph with a laminar family C = {S'j}^]^ of vertex-sets having bounds ^; and the goal is to 
compute a spanning tree of minimum cost that contains at most 6/ edges from 5{Si) for each i G [m]. 

The motivation behind this problem is in designing a network where there is a hierarchy (i.e. laminar 
family) of service providers that control nodes (i.e. vertices). The number of edges crossing the boundary 
of any service provider (i.e. its vertex-cut) represents some cost to this provider, and is therefore limited. 
The laminar MCST problem precisely models the question of connecting all nodes in the network while 
satisfying bounds imposed by all the service providers. 

From a theoretical viewpoint, cut systems induced by laminar families are well studied, and are 
known to display rich structure. For example, one-way cut-incidence matrices are matrices whose 
rows are incidence vectors of directed cuts induced by the vertex-sets of a laminar family; It is well 
known (e.g., see [j53l) that such matrices are totally unimodular. Using the laminar structure of degree- 
constraints and the iterative relaxation framework, we obtain the following main result, and present its 
proof in Section [2l 

Theorem 1 There is a polynomial time (1, 6 + 0(log n)) bicriteria approximation algorithm for laminar 
MCST. That is, the cost is no more than the optimum cost and the degree violation is at most additive 
0(log n). This guarantee is relative to the natural LP relaxation. 

This guarantee is substantially stronger than what follows from known results for the general min- 
imum crossing spanning tree (MCST) problem: where the degree bounds could be on arbitrary edge- 
subsets El, . . . , Em- In particular, for general MCST a (1, 6 + A — 1) H [191 is known where A is 
the maximum number of degree-bounds an edge appears in. However, this guarantee is not useful for 
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laminar MCST as A can be as large as in this case. If a multiplicative factor in the degree violation 
is allowed, Chekuri et al. fTTl recently gave a very elegant (l,(l + e)6 + 0(^ log m)) guarantee (which 
subsumes the previous best (O(logn), O(logm) 6) 161 result). However, these results also cannot be 
used to obtain a small additive violation, especially if b is large. In particular, both the results 161 [HI for 
general MCST are based on the natural LP relaxation, for which there is an integrality gap ofb + Q,{^/n) 
even without regard to costs and when m = 0{n) ll30l (see also Section ll!2l ). On the other hand. The- 
orem [T] shows that a purely additive O(logn) guarantee on degree (relative to the LP relaxation and 
even in presence of costs) is indeed achievable for MCST, when the degree-bounds arise from a laminar 
cut-family. 

The algorithm in Theorem[T]is based on iterative relaxation and uses two main new ideas. Firstly, we 
drop a carefully chosen constant fraction of degree-constraints in each iteration. This is crucial as it can 
be shown that dropping one constraint at a time as in the usual applications of iterative relaxation can 
indeed lead to a degree violation of ri(A). Secondly, the algorithm does not just drop degree constraints, 
but in some iterations it also generates new degree constraints, by merging existing degree constraints. 

All previous applications of iterative relaxation to constrained network design treat connectivity and 
degree constraints rather asymmetrically. While the structure of the connectivity constraints of the under- 
lying LP is used crucially (e.g., in the ubiquitous uncrossing argument), the handling of degree constraints 
is remarkably simple. Constraints are dropped one by one, and the final performance of the algorithm is 
good only if the number of side constraints is small (e.g., in recent work by Grandoni et al. fT6l ). or if 
their structure is simple (e.g., if the 'frequency' of each element is small). In contrast, our algorithm for 
laminar- MCST exploits the structure of degree constraints in a non-trivial manner. 

Hardness Results. We obtain the following hardness of approximation for the general MCST problem 
(and its matroid counterpart). In particular this rules out any algorithm for MCST that has additive 
constant degree violation, even without regai^d to costs. 

Theorem 2 Unless MV has quasi-polynomial time algorithms, the MCST problem admits no polyno- 
mial time 0(log" m) additive approximation for the degree bounds for some constant a > 0; this holds 
even when there are no costs. 

The proof for this theorem is given in Section|3l and uses a a two-step reduction from the well-known 
Label Cover problem. First, we show hardness for a uniform matroid instance. In a second step, we then 
demonstrate how this implies the result for MCST claimed in Theorem |2] 

Note that our hardness bound nearly matches the result obtained by Chekuri et al. in ITTI . We note 
however that in terms of purely additive degree guarantees, a large gap remains. As noted above, there is 
a much stronger lower bound of 6 + O(y^) for LP-based algorithms l30l (even without regard to costs), 
which is based on discrepancy. In light of the small number of known hardness results for discrepancy 
type problems, it is unclear how our bounds for MCST could be strengthened. 

An interesting consequence of the hai^dness result in Theorem |2] is for the robust (or min-max) k- 
median problem HI. In this problem, there are m different client-sets in a metric and the goal is to 
open k facilities that are simultaneously good (in terms of the A;-median objective) for all the client-sets. 
Anthony et al. HI obtained a logarithmic approximation algorithm for this problem, and showed that it 
is hard to approximate better than factor 2. The following result shows that the robust fc-median problem 
is indeed harder to approximate than usual /c-median, for which 0(l)-approximations are known ITIO. 
We present its proof in Section [TTl 

Corollary 3 Robust k-median is 0(log° m)-hard to approximate even on uniform metrics (for some 
fixed constant a> Q), assuming MV does not have quasi-polynomial time algorithms. 
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Degree Bounds in More General Settings. We consider crossing versions of other classic combina- 
torial optimization problems, namely contra-polymatroid intersection and lattice polyhedra |[29l . 

Definition 4 (Minimum crossing contra-polymatroid intersection problem) Let ri,r2 2^ ^ Z be 
two supermodular functions, c : E ^ ]R_|_ and {Ei}i^£ be a collection of subsets ofE with corresponding 
bounds {bi}i£i. Then the goal is to minimize: 

{c^x I x{S) > max{ri{S),r2{S)},\/ S C E; 
x{Ei) <h, Mi^I] x£ {0,1}^}. 

In particular, this definition captures the degree-bounded version of spanning-set intersection in two 
matroids (for eg. the bipartite edge-cover problem). We note that this definition does not capture alternate 
notions of matroid intersection, such as intersection of bases in two matroids; hence it does not apply to 
the degree-bounded arborescence problem. 

Let A = maXggE \ {i G [m] \ e G Ei}\ be the largest number of sets Ei that any element of E 
belongs to, and refer to it as frequency. The proof of this theorem can be found in Section|4] 

Theorem 5 Any optimal basic solution x* of the linear relaxation of the minimum crossing contra- 
polymatroid intersection problem can be rounded into an integral solution x such that: 

x{S) > max{ri(S),r2(S)}, VS C E; x{Ei) < 26^ + A - 1, Vi G /; and c^x < 2c^x*. 

The algorithm for this theorem again uses iterative relaxation, and its proof is based on a 'fractional 
token' counting argument similar to the one used in 141. We also observe that the natural iterative relax- 
ation steps are insufficient to obtain a better approximation guarantee. 

Crossing Lattice Polyhedra. Classical lattice polyhedra form a unified framework for various discrete 
optimization problems and go back to Hoffman and Schwartz [|171 who proved their integrality. They are 
polyhedra of type 

{x G [0, 1]-^ I x{p{S)) > r{S), yS G T} 

where is a consecutive submodular lattice, p : — 2^ is a mapping from T to subsets of the ground- 
set E, and r G M.-^ is supermodular. A key property of lattice polyhedra is that the uncrossing technique 
can be applied which turns out to be crucial in almost all iterative relaxation approaches for optimization 
problems with degree bounds. We refer the reader to |[29l for a more comprehensive treatment of this 
subject. 

We generalize our work further to crossing lattice polyhedra which arise from classical lattice polyhe- 
dra by adding "degree-constraints" of the form Cj < x{Ei) < bi for a given collection {Ei ^ E' | i g /} 
and lower and upper bounds a,b & M^. We mention two (of several) examples which are covered by 
this model: 

Example 1: Crossing matroid basis. Here T = 2^, p is the identity map, and the partial order in 
is the canonical one that is induced by set inclusion. Function r : 2^ — )• N is defined as r{S) = 
rank(V) — rank(V \ S); where E is the ground-set of the matroid and rank is its rank function. The 
crossing matroid basis problem finds the minimum cost basis in the matroid satisfying degree bounds. 

Example 2: Crossing planar min cut. Let G = {V, E) be a (directed or undirected) s, t-planar graph 
(along with an embedding) with s, i G Here elements of F correspond to s-t paths in G (p maps each 

'in an earlier version of the paper ||5j, we had incorrectly claimed that our result extends to degree-bounded arborescence. 
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element of T to the edge-set of that s — t path), and the partial order in T relates paths where one is 
below/above the other in the planar embedding of G. The rank function is the constant all-ones function. 
The crossing planar min-cut problem involves finding a minimum cost s — t cut in G that obeys the 
degree bounds. 

We can show that the standard LP relaxation for the general crossing lattice polyhedron problem is 
weak; in Section lSTTl we give instances of crossing planar min-cut (i.e., Example 2 above) where the LP- 
relaxation is feasible, but any integral solution violates some degree-bound by VL[^/n). For this reason, 
we henceforth focus on a restricted class of crossing lattice polyhedra in which the underlying lattice 
(-^, <) satisfies the following monotonicity property 

(*) S<T ^ < |p(T)| V5,TG^. 

We obtain the following theorem whose proof is given in Section |5] 

Theorem 6 For any instance of the crossing lattice polyhedron problem in which T satisfies property 
(*), there exists an algorithm that computes an integral solution of cost at most the optimal, where all 
rank constraints are satisfied, and each degree bound is violated by at most an additive 2A — 1. 

We note that the above property (*) is satisfied for matroids, and hence Theorem |6] matches the 
previously best-known bound |[T9l for degree bounded matroids (with both upper/lower bounds). Also 
note that property (*) holds whenever F is ordered by inclusion. In this special case, we can improve 
the result to an additive A — 1 approximation if only upper bounds are given. 

1.2 Related Work 

As mentioned earlier, the basic bounded-degree MST problem has been extensively studied E] HI [151 
[21] |22j |27l HH |3ll. The iterative relaxation technique for degree-constrained problems was developed 
in El HI. 

MCST was first introduced by Bilo et al. O, who presented a randomized-rounding algorithm that 
computes a tree of cost 0(log n) times the optimum where each degree constraint is violated by a mul- 
tiplicative O(logn) factor and an additive O(logm) term. Subsequently, Bansal et al. H gave an algo- 
rithm that attains an optimal cost guarantee and an additive A — 1 guai^antee on degree; recall that A 
is the maximum number of degree constraints that an edge lies in. This algorithm used iterative relax- 
ation as its main tool. Recently, Chekuri et al. [11] obtained an improved (l, (1 + e)b + 0{j logm)) 
approximation algorithm for MCST, for any e > 0; this algorithm is based on pipage rounding. 

The minimum crossing matroid basis problem was introduced in |[T9l . where the authors used itera- 
tive relaxation to obtain (1)(1,6 + A — 1) -approximation when there are only upper bounds on degree, 
and (2) (1, b + 2A — 1) -approximation in the presence of both upper and lowed degree-bounds. The ifTTI 
result also holds in this matroid setting. [191 also considered a degree-bounded version of the submodular 
flow problem and gave a (1, 6 + 1) approximation guarantee. 

The bounded-degree arborescence problem was considered in Lau et al. ll24l . where a (2, 26 + 2) 
approximation guarantee was obtained. Subsequently Bansal et al. H| designed an algorithm that for 
any < e < 1/2, achieves a (1/e, — e) + 4) approximation guarantee. They also showed that 
this guarantee is the best one can hope for via the natural LP relaxation (for every < e < 1/2). In 
the absence of edge-costs, |4j gave an algorithm that violates degree bounds by at most an additive two. 
Recently Nutov [261 studied the arborescence problem under weighted degree constraints, and gave a 
(2, 56) approximation for it. 
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Lattice polyhedra were first investigated by Hoffman and Scliwartz lITTl and the natural LP relaxation 
was shown to be totally dual integral. Even though greedy-type algorithms are known for all examples 
mentioned earlier, so far no combinatorial algorithm has been found for lattice polyhedra in general. 
Two-phase greedy algorithms have been established only in cases where an underlying rank function 
satisfies a monotonicity property lfT4l . llT2l . 



2 Crossing Spanning Tree with Laminar degree bounds 

In this section we prove Theorem [T] by presenting an iterative relaxation-based algorithm with the stated 
performance guarantee. During its execution, the algorithm selects and deletes edges, and it modifies 
the given laminar family of degree bounds. A generic iteration starts with a subset F of edges already 
picked in the solution, a subset E of undecided edges, i.e., the edges not yet picked or dropped from the 
solution, a laminai^ family C on V, and residual degree bounds b{S) for each 5 € £. 

The laminar family C has a natural forest-like structure with nodes corresponding to each element of 
C. A node 5 G £ is called the parent of node C G £ if 5 is the inclusion-wise minimal set in £ \ {C} 
that contains C; and C is called a child of S. Node € £ is called a grandchild of node S € £ if S" is 
the parent of D's parent. Nodes S,T ^ C are siblings if they have the same parent node. A node that 
has no parent is called root. The level of any node 5 G £ is the length of the path in this forest from S to 
the root of its tree. We also maintain a linear ordering of the children of each £-node. A subset B C C 
is called consecutive if all nodes in B are siblings (with parent S) and they appear consecutively in the 
ordering of 5's children. In any iteration {F, E, C,b), the algorithm solves the following LP relaxation 
of the residual problem. 



min ^ CeXe (1) 

eeE 

S.t. x{E{V)) = \V\ - |F| - 1 

x{E{U)) <\U\-\F{U)\~1 yUdV 
x{5e{S)) < b{S) ySeC 
Xe>0 \fee E 



For any vertex-subset W <^ V and edge-set H, we let H{W) := {{u,v) H \ u,v ^ W} denote 
the edges induced on W; and 6h{W) := {{u,v) H \ u ^ W, v ^ W} the set of edges crossing W. 
The first two sets of constraints are spanning tree constraints while the third set corresponds to the degree 
bounds. Let x denote an optimal extreme point solution to this LP. By reducing degree bounds h{S), if 
needed, we assume that x satisfies all degree bounds at equality (the degree bounds may therefore be 
fractional-valued). Let a := 24. 

Definition 7 An edge e £ E is said to be local for S G C if e has at least one end-point in S but is 
neither in E{C) nor in S{C) Pi 5{S)for any grandchild C of S. Let local(S') denote the set of local edges 
for S. A node S G C is said to be good //'|local(S')| < a. 

Figure [T] shows a set S, its children Bi and B2, and grand-children Ci, . . . ,0^; edges in local(5) 
aie drawn solid, non-local ones are shown dashed. 

Initially, E is the set of edges in the given graph, F 0, £ is the original laminar family of vertex 
sets for which there are degree bounds, and an arbitrary linear ordering is chosen on the children of each 
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Figure 1: Example of local edges. 



node in £. In a generic iteration {F, E, C,b), the algorithm performs one of the following steps (see also 
Figure 121): 

1. lfxe = l for some edge e e E then F ^ FU {e}, E ^ E\ {e}, and set 6(5') ^ b{S) - 1 for 
all 5 G £ with e G 6{S). 

2. If Xe = for some edge e G E then E ^ E \ {e}. 

3. DropN: Suppose there at least \C\/A good non-leaf nodes in C Then either odd-levels or even- 
levels contain a set 7W C £ of |£|/8 good non-leaf nodes. Drop the degree bounds of all children 
of A4 and modify C accordingly. The ordering of siblings also extends naturally. 

4. DropL: Suppose there are more than |/2|/4 good leaf nodes in C, denoted by M. Then partition J\f 
into parts coiTcsponding to siblings in C. For any part {A'^i, • • • , Nk} C J\f consisting of ordered 
(not necessarily contiguous) children of some node 5: 

(a) Define Mi = iV2i-i U N2i for all 1 < i < [k/2\ (if k is odd is not used). 

(b) Modify C by removing leaves {A'^i, • • • , iV^} and adding new leaf-nodes {Mi, • • • , M|^fc/2J } 
as children of 5 (if k is odd Nk is removed). The children of S in the new laminar family 
ai^e ordered as follows: each node Mj takes the position of either N2i~i or N2i, and other 
children of S are unaffected. 

(c) Set the degree bound of each Mi to h{Mi) = b{N2i^i) + b{N2i). 

Assuming that one of the above steps applies at each iteration, the algorithm terminates when E = % 
and outputs the final set F as a solution. It is clear that the algorithm outputs a spanning tree of G. An 
inductive argument (see e.g. |[24l ) can be used to show that the LP ([T]) is feasible at each each iteration 
and c{F) + Zcur < where Zq is the original LP value, Zcur is the current LP value, and F is the chosen 
edge-set at the current iteration. Thus the cost of the final solution is at most the initial LP optimum Zq- 
Next we show that one of the four iterative steps always applies. 

Lemma 8 In each iteration, one of the four steps above applies. 

Proof: Let x* be the optimal basic solution of ([U, and suppose that the first two steps do not apply. 
Hence, we have < x* < 1 for all e ^ E. The fact that x* is a basic solution together with a standai^d 
uncrossing argument (e.g., see implies that x* is uniquely defined by 

x{E{U)) = \U\ - \F{U)\ - 1 VC/GcS, and x{6EiS)) = b{S), ySeC, 
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Figure 2: Examples of the degree constraint modifications DropN and DropL. 



where 5 is a laminar subset of the tight spanning tree constraints, and C is a subset of tight degree 
constraints, and where \E\ = \S\ + \C'\. 

A simple counting argument (see, e.g., 11311 ) shows that there are at least 2 edges induced on each 
S £ S that are not induced on any of its children; so 2\S\ < \E\. Thus we obtain \E\ < 2\C'\ < 2\C\. 

From the definition of local edges, we get that any edge e = (n, v) is local to at most the following 
six sets: the smallest set Si £ C containing u, the smallest set S2 ^ C containing v, the parents Pi 
and P2 of Si and S2 resp., the least-common-ancestor L of Pi and P2, and the parent of L. Thus 
^^g^ |local(S')| < Q\E\. From the above, we conclude that I local(S') I < 12|£|. Thus at least 

\C\/2 sets S e C must have |local(S')| < a = 24, i.e., must be good. Now either at least |>C|/4 of them 
must be non-leaves or at least \C\/A of them must be leaves. In the first case, step 3 holds and in the 
second case, step 4 holds. ■ 

It remains to bound the violation in the degree constraints, which turns out to be rather challenging. 
We note that this is unlike usual applications of iterative rounding/relaxation, where the harder part is in 
showing that one of the iterative steps applies. 

It is clear that the algorithm reduces the size of C by at least |£|/8 in each DropN or DropL iteration. 
Since the initial number of degree constraints is at most 2n — 1, we get the following lemma. 

Lemma 9 The number of drop iterations (DropN and DropL) is T := 0(log n). 

Performance guarantee for degree constraints. We begin with some notation. The iterations of the 
algorithm are broken into periods between successive drop iterations: there are exactly T drop-iterations 
(Lemma|9ll. In what follows, the t-th drop iteration is called round t. The time t refers to the instant just 
after round t; time refers to the start of the algorithm. At any time t, consider the following parameters. 

• Ct denotes the laminar family of degree constraints. 

• Et denotes the undecided edge set, i.e., support of the current LP optimal solution. 

• For any set B of consecutive siblings in Ct, Bnd(i3,t) = Y^NeB^^^) equals the sum of the 
residual degree bounds on nodes of B. 

• For any set B of consecutive siblings in Ct, lnc(fi, t) equals the number of edges from (U ngB^) 
included in the final solution. 
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Recall that b denotes the residual degree bounds at any point in the algorithm. The following lemma 
is the main ingredient in bounding the degree violation. 



Lemma 10 For any set B of consecutive siblings in Ct (at any time t), Inc(i3,t) < Bnd(B, t) + 4a • (T — 
t). 

Observe that this implies the desired bound on each original degree constraint S: using t = and 
B = {S}, the violation is bounded by an additive 4a • T term. 

Proof: The proof of this lemma is by induction onT — t. The base case t = T is trivial since the only 
iterations after this correspond to including 1 -edges: hence there is no violation in any degree bound, 
i.e. \nc{{N},T) < b{N) for all N G Ct- Hence for any B Q C, \nc{B,T) < Y.Net3 lnc({iV},T) < 
j:NeBKN) = Bnd{B,T). 

Now suppose t < T, and assume the lemma for i + 1. Fix a consecutive B Ct- We consider 
different cases depending on what kind of drop occurs in round t + 1. 

DropN round. Here either all nodes in B get dropped or none gets dropped. 

Case 1: None of B is dropped. Then observe that B is consecutive in as well; so the inductive 
hypothesis implies Inc(i3, < Bnd(i3, t+l)+4a-(T— t— 1). Since the only iterations between round 
t and round i + 1 involve edge-fixing, we have \nc{B, t) < Bnd{B, t) — Bnd(;B, t + l) + \nc{B, t + 1) < 
Bnd{B, t) + 4a-{T-t-l)< Bnd(B, t) + 4a ■ (T - t). 

Case 2: All ofB is dropped. Let C denote the set of all children (in Ct) of nodes in B. Note that C 
consists of consecutive siblings in and inductively lnc(C,t + l) < Bnd(C, t + 1) + 4a • (T — f — 1). 
Let S ^ Ct denote the parent of the ;B-nodes; so C are grand-children of 5 in Ct. Let x denote the 
optimal LP solution just before round t+l (when the degree bounds are still given by Ct), and H = 
Et+i the support edges of x. At that point, we have b{N) = x{6{N)) for all N £ B U C. Also let 
Bnd'(B,t + 1) := '^N^s^i-^) bounds on i3-nodes just before round t + l. Since S 

is a good node in round t + 1, |Bnd'(^,t + 1) - Bnd(C,i + 1)| = | YIngbK^) " EMeC^(^^)l = 

1 ^AfeZ? ^i^i^)) ~ J2Aiec — The last inequality follows since S is good; the factor of 

2 appeai^s since some edges, e.g., the edges between two children or two grandchildren of S, may get 
counted twice. Note also that the symmetric difference of dni'^NeB^) and dni^AiecM) is contained 
in local(S'). Thus dni^NeB^) and Sni^MecM) differ in at most a edges. 

Again since all iterations between time t and i + 1 are edge-fixing: 

\nc{B,t) < Bnd{B,t) -Bnd'iB, t + l) + \6H{UNeBN)\5H{UM^cM)\ 
+lnc(C,t + 1) 

< Bnd{B, t) - Bnd'(S, t + l) + a + lnc(C, t + 1) 

< Bnd(i3, t) - Bnd'{B, t + l) + a + Bnd(C, t + 1) + 4a • (T - t - 1) 

< Bnd(i3, t) - Bnd'{B, t + l) + a + Bnd\B, t + 1) + 2a + 4a ■ [T - t - 1) 

< Bnd{B,t) + 4a-{T-t) 

The first inequality above follows from simple counting; the second follows since dni^NeB^) 
and dni^MecM) differ in at most a edges; the third is the induction hypothesis, and the fourth is 
Bnd(C, t + l) < Bnd'(i3, t + 1) + 2a (as shown above). 

DropL round. In this case, let S be the parent of S-nodes in Ct, and = {A'^i, • • • , Np} be all the 
ordered children of S, of which S is a subsequence (since it is consecutive). Suppose indices 1 < 7r(l) < 
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7r(2) < • • • < TT{k) < p correspond to good leaf-nodes in M. Then for each 1 < i < [k/2\, nodes 
-^7r(2i-i) ^rid N^(^2i) merged in this round. Let {vr(i) | e < i < /} (possibly empty) denote the 
indices of good leaf-nodes in B. Then it is clear that the only nodes of B that may be merged with nodes 
outside B are N^(^f.^ and iV^(/) ; other ;B-nodes are either not merged or merged with another S-node. 
Let C be the inclusion-wise minimal set of children of S in Ct+i s.t. 

• C is consecutive in 

• C contains all nodes of B\ {iV^(i)}f=i, and 

• C contains all new leaf nodes resulting from merging two good leaf nodes of B. 

Note that Um^^qM consists of some subset of B and at most two good leaf-nodes inM \ B. These 
two extra nodes (if any) are those merged with the good leaf-nodes A^,r(e) and -/V^(/) of B. Again let 
Bnd'(;B,t + 1) := J^NeB^^) denote the sum of bounds on B just before drop round t + 1, when 
degree constraints are Ct- Let H = Et+i be the undecided edges in round t + 1. By the definition of 
bounds on merged leaves, we have Bnd(C, t + l)< Bnd'(5, t + l) + 2a. The term 2a is present due to 
the two extra good leaf-nodes described above. 

Claim 11 Wehave \5H{^NeBN) \ 5h{^M(lcM)\ < 2a. 

Proof: We say that N ^ J\f is represented in C if either G C or iV is contained in some node of C. 
Let V be set of nodes of B that are not represented in C and the nodes of J\f\B that ai^e represented in C. 
Observe that by definition of C, the set V C {N^(_e-i)i ^^{e)^ -^7r(/); -^7r(/+i)}; iri fact it can be easily 
seen that \V\ < 2. Moreover V consists of only good leaf nodes. Thus, we have | ULev Sh{L)\ < 2a. 
Now note that the edges in Sni^NeBN) \ 6h{^m&cM) must be in LlLev^HiL). This completes the 
proof. ■ 

As in the previous case, we have: 

\nc{B,t) < Bnd{B,t) -Bnd'{B,t + l) + \SH{UNeBN)\5H{UMecM)\ 
+lnc(C,t + 1) 

< Bnd(e, t) - Bnd'{B, t + 1) + 2a + lnc(C, t + 1) 

< Bnd(i3, t) - Bnd'(i3, t + 1) + 2a + Bnd(C, t + 1) + 4a • (T - t - 1) 

< Bnd{B, t) - Bnd'{B, t + 1) + 2a + Bnd'{B, t + 1) + 2a + 4a ■ {T - t - 1) 
= Bnd{B,t) + 4a-{T-t) 

The first inequality follows from simple counting; the second uses Claim[TTl the third is the induction 
hypothesis (since C is consecutive), and the fourth is Bnd(C, t + 1) < Bnd'{B, t + 1) + 2a (from above). 
This completes the proof of the inductive step and hence Lemma [TOl ■ 



3 Hardness Results 

In this section we prove Theorem |2l i.e. unless J\fV has quasi-polynomial time algorithms, there is 
no polynomial time O(log'^m) additive approximation for degree bounds for the minimum crossing 
spanning tree problem, where c > is some universal constant. This result also holds in the absence 
of edge-costs. We note that this hardness result only holds for the general MCST problem, and not the 
laminar MCST addressed earlier. The first step to proving this result is a hardness for the more general 
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minimum crossing matroid basis problem: given a matroid on a ground set V of elements, a cost 
function c : V ^ ]R+, and degree bounds specified by pairs {{Ei,bi)}^^ (where each Ei V and 
bi G N), find a minimum cost basis I in A4 such that | / n < hi for all i G [m]. 

Theorem 12 Unless MV has quasi-polynomial time algorithms, the unweighted minimum crossing ma- 
troid basis problem admits no polynomial time O(log^ m) additive approximation for the degree bounds 
for some fixed constant c > 0. 

Proof: We reduce from the label cover problem [2|. The input is a graph G = {U, E) where the vertex 
set U is partitioned into pieces C/i , • • • , each having size q, and all edges in E are between distinct 
pieces. We say that there is a superedge between Ui and Uj if there is an edge connecting some vertex in 
Ui to some vertex in Uj. Let t denote the total number of superedges; i.e., 



The goal is to pick one vertex from each part {Ui]^^^ so as to maximize the number of induced 
edges. This is called the value of the label cover instance and is at most t. 

It is well known that there exists a universal constant 7 > 1 such that for every k G N, there is a 
reduction from any instance of SAT (having size N) to a label cover instance {G = ([/, E), q, t) such 
that: 

• If the SAT instance is satisfiable, the label cover instance has optimal value t. 

• If the SAT instance is not satisfiable, the label cover instance has optimal value < t/j'^. 

• |G| = iY^C^), q = 2^, \E\ < t^, and the reduction runs in time N'^C'). 

We consider a uniform matroid A4 with rank t on ground set E (recall that any subset of t edges is 
a basis in a uniform matroid). We now construct a crossing matroid basis instance Z on A4. There is a 
set of degree bounds corresponding to each i G [n] : for every collection C of edges incident to vertices 
in Ui such that no two edges in G are incident to the same vertex in Ui, there is a degree bound in X 
requiring at most one element to be chosen from C. Note that the number of degree bounds m is at most 
\E\i < ArO(fc2''). ji^g following claim links the SAT and crossing matroid instances. 

Claim 13 [ Yes instance ] If the SAT instance is satisfiable, there is a basis (i.e. subset B C E with 
\B\ = t) satisfying all degree bounds. 

[ No instance ] If the SAT instance is unsatisfiable, every subset B' Q E with \B'\ > t/2 violates some 
degree bound by an additive p = 7'^/^/\/2. 

Proof: Observe that if the original SAT instance is satisfiable, then the matroid A4 contains a basis 
obeying all the degree bounds: namely the t edges T* (1 E covered in the optimal solution to the label 
cover instance. This is because if we consider any Ui, then all the T*-edges having a vertex in Ui as their 
endpoint, have the same endpoint. Thus, for any degree bound corresponding to collection G (as defined 
above), at most one T*-edge can lie in G. 

Now consider the case that the SAT instance is unsatisfiable. Let B' (1 E he. any subset with 
\B'\ > t/2. We claim that B' contains at least p = 7'^/^/\/2 edges from some degree-constrained 
set of edges. Suppose (for a contradiction) that \B' Ci G\ < p for each degree constraint G. This 
means that each part {Ui}f^i contains fewer than p vertices that ai^e incident to edges B. For each part 
i G [n], let Wi C Ui denote the vertices incident to edges of B; note that \ Wi\ < p. Consider the label 
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cover solution obtained as follows. For each i G [n], choose one vertex from Wi independently and 
uniformly at random. Clearly, the expected number of edges in the resulting induced subgraph is at least 
iB'l/p"^ > 2^ = t/^''. This contradicts the fact that the value of label cover instance is strictly less than 
t/^f''. m 

The steps described in the above reduction can be done in time polynomial in m and Also, 
instead of randomly choosing vertices from the sets Wi, we can use conditional expectations to derive 
a deterministic algorithm that recovers at least t/p^ edges. Setting k = B(loglog A^) (recall that is 
the size of the original SAT instance), we obtain an instance of bounded-degree matroid basis of size 
max{m, |G|} = N^"^"" ^ and p = log^ N, where a,b > are constants. Note that log m = log""*'^ N, 
which implies p = log'^ m for c = > 0, a constant. Thus it follows that for this constant c > the 
bounded-degree matroid basis problem has no polynomial time 0(log'^?Ti) additive approximation for 
the degree bounds, unless MV has quasi-polynomial time algorithms. ■ 

We now prove Theorem |2l 
Proof: [Proof of Theorem |2l We show how the bases of a uniform matroid can be represented in a 
suitable instance of the crossing spanning tree problem. Let the uniform matroid from Theorem [12] 
consist of e elements and have rank t < e; recall that t > ^/e and clearly m < 2^. We construct a graph 
as in Figure [3l with vertices ui , • • • ,Ve con^esponding to elements in the uniform matroid. Each vertex 
Vi is connected to the root r by two vertex-disjoint paths: {vi,Ui,r) and {vi,Wi,r). There are no costs in 
this instance. Corresponding to each degree bound (in the uniform matroid) of b{C) on a subset C C [e], 
there is a constraint to pick at most |C| + b{C) edges from 6{{ui \ i G C}). Additionally, there is a 
special degree bound of 2e — t on the edge-set E' = [Jl^i S{wi); this corresponds to picking a basis in 
the uniform matroid. 




Vi 



Figure 3: The crossing spanning tree instance used in the reduction. 

Observe that for each i £ [e], any spanning tree must choose exactly three edges amongst {(r, Ui), {ui,Vi), 
(r, Wi), {wi,Vi)}, in fact any three edges suffice. Hence every spanning tree T in this graph corresponds 
to a subset X C [e] such that: (I) T contains both edges in 6{ui) and one edge from 5{wi), for each 
i £ X, and (II) T contains both edges in 6{wi) and one edge from 6{ui) for each i G [e] \ X. 

From Theorem [T2I for the crossing matroid problem, we obtain the two cases: 

Yes instance. There is a basis B* (i.e. B* (Z [e], \B*\ = t) satisfying all degree bounds. Consider the 
spanning tree 

T* = {(r, Ui), {ui,Vi), (r, Wi) \ i € B*} [J{{r, Wi), {ui,Wi), (r, Ui) \ i £ [e] \ B*}. 

Since B* satisfies its degree-bounds, T* satisfies all degree bounds derived from the crossing matroid 
instance. For the special degree bound on E', note that \T* n E'\ = 2e — \B*\ = 2e — t; so this is also 
satisfied. Thus there is a spanning tree satisfying all the degree bounds. 
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No instance. Every subset B' C [e] with \B'\ > t/2 (i.e. neai^ basis) violates some degree bound by 
an additive p = Q{log'^ m) term, where c > is a fixed constant. Consider any spanning tree T that 
corresponds to subset X C [e] as described above. 

1. Suppose that \X\ < t/2; then we have |T n E'\ = 2e — \X\ > 2e — t + |, i.e. the special degree 
bound is violated by t/2 > Q{y/e) = r2(log^/^ m). 

2. Now suppose that \X\ > t/2. Then by the guarantee on the no-instance, T violates some degree- 
bound derived from the crossing matroid instance by additive p. 

Thus in either case, every spanning tree violates some degree bound by additive p = r2(log'^ m). 

By Theorem [121 it is hard to distinguish the above cases and we obtain the corresponding hai^dness 
result for crossing spanning tree, as claimed in Theorem |2] H 

3.1 Hardness for Robust fc-median 

Another interesting consequence of Theorem [12] is for the robust A;-median problem |[T1. Here we are 
given a metric {V,d), m client-sets {Si C F}™^, and bound k; the goal is to find a set F C y of A; 
facilities such that the worst-case connection cost (over all client-sets) is minimized, i.e. 

mill nmx div, F). 
FCV,\F\=k 1=1 ^ ^ 

Above d{v, F) denotes the shortest distance from v to any vertex in F. Anthony et al. m gave an 
0(log m+log /c)-approximation algorithm for robust /c-median, and showed that it is hard to approximate 
better than factor two. At first sight this problem may seem unrelated to crossing matroid basis. However 
using Theorem [T2j we obtain the poly-logarithmic hardness result stated in Corollary [3] 
Proof: Recall that in a uniform metric, the distance between every pair of vertices is one. In this case 
the robust A;-median problem can be rephrased as: 

min rnax IS'j \ F|, where jS'j C T/)^, are the client-sets. 

FCV,\F\=k i=l ' * ^ " ^ ' - ^'-^ 

The hard instances of crossing matroid basis in Theorem [12] are in fact for uniform matroids where 
every degree upper-bound equals one. i.e. there is a ground-set V, degree bounds given by {Ei C V}^^, 
and rank t; the goal is to find (if possible) a subset / C 1/ with |/| = t such that |/P| < 1 for all 
i G [m]. Theorem [T2] showed that it is hard to distinguish the following cases: (Yes-case) there is 
some / C y with |/| = t and maxjg[f„] \I n Ei\ < 1; and (No-case) for every I <^ V with |/| = t, 
maxjg[m] \I nEi\> p := ^{log^'m). 

These hard instances naturally correspond to the robust /c-median problem on uniform metric V, 
client-sets {Ei C y}™^, and bound k = \V\ — t. It is clear that the robust /c-median objective is at 
most one in the Yes-case, and at least p in the No-case. Thus we obtain a multiplicative p hardness of 
approximation for robust A;-median on uniform metrics. This proves Corollary [3] ■ 

3.2 Integrality Gap for general MCST 

We now present the h + Q.{y/ri) integrality gap instance for minimum crossing spanning tree. While 
such gaps instances ai^e easy to obtain if one allows m to be super-polynomially lai^ge (for example, by 
setting a degree bound for each subset of edges), the nice property of the example here is that m is quite 
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small, in fact m = 0{n). This result is due to Mohit Singh ll30l . we thank him for letting us present the 
example here. 

The graph is the same as the one used for the hardness result. The vertex-set is {r} [j{vi,Ui, Wi}|=i 
so n = 3e + 1. The edges are {(r, Ui) \ i e [e]} U {{vi,Ui) \ i G [e]} and {(r, Wi) \ i G [e]} U {{vi,Wi) \ 
i G [e]}. See also Figure[3] There are no costs in this instance. 

The 'degree bounds' for the MCST instance are derived from the lower bound for the discrepancy 
problem lITOl . From discrepancy theory there exists a collection {Sj C [e]}^^^ of subsets such that, 



Above X = [e] \ X as usual, and p = Q.{y/e) = Q.{y/n). In other words, for every way of partitioning 
[e], there is some set Sj such that the partition induced on Sj has a large imbalance. There are m = 2e 
degree bounds, defined as follows. For each j G [e] there is a bound of \Sj \ + [|5j|/2] on each of the 
edge-sets Uj = Ui^Sj^Ui) = {(r,itj), {ui,Vi)}i(zs^, and Wj = Ui(zs^6{wi) = {{r,Wi), {wi,Vi)}i^Sj- 

Consider the fractional solution to the natural LP relaxation that sets each edge to value 3/4. It is 
easily seen that it is indeed a fractional spanning tree and satisfies all the degree bounds. 

On the other hand, we claim that any integer solution must violate some degree bound by additive 
^ — 1. Note that every spanning tree T in this graph corresponds to a subset X C [e] such that: (I) T 
contains both edges in 6{ui) and one edge from 6{wi), for each i £ X, and (II) T contains both edges in 
6{iUi) and one edge from 6{ui) for each i G X. The number of edges used by tree T in the degree-bounds 
(for each j G [e]) are: 

• \Tr\Uj\=2\XnSj\ + \Xn Sj\ = \Sj\ + \XnSj\, and 

• \TnWj\ = \x nSj\ + 2\xnSj\ = \Sj\ + \xnSj\. 

From the discrepancy instance, it follows that max^^^^ | |X n S'jl — |X n S'jl | > p; let A; be the index 
achieving this maximum. Then we have: 



4 Minimum Crossing Contra-Polymatroid Intersection 

In this section we consider the crossing contra-polymatroid intersection problem (see Definition |4| and 
prove Theorem [5] The algorithm (given as Algorithm [B for this problem is based on iteratively relaxing 
the following natural LP relaxation. 



max \ \X r\ Sj\ - \X r\ Sj\ \ > p, for every X C [e]. 



max{|rn[/,.|, \TnWk\} = \Sk\+meix{\XnSk\, \xnSk\}>\Sk\ + ^-^ + ^. 



Thus the degree-bound for either Uk or is violated by additive f — 1. 




x{SnE') > ri(5) - |FnS| 
x{SnE') > r2{S) -\FnS\ 
x{E, n E') < b'i 
< Xe < 1 



yi£W 

Ve G E'. 



yscE 



ysoE 
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At a generic iteration, E' (1 E denotes the set of unfixed elements, F Q E the set of chosen elements 
(recall that E denotes the groundset of the instance), W <^ I the set of remaining degree bounds, and 
6^ (for each i € W) the residual degree-bound in the i^^ constraint. Observe that this LP can indeed 
be solved in polynomial time by the Ellipsoid algorithm: the separation oracle for the first two sets of 
constraints involve submodular function minimization for the two functions gi{S) = x{S n E') + \ S Ci 
F\ — ri{S) (with i = 1,2). The resulting fractional solution can then be converted to an extreme point 
solution of no larger cost, as described in Jain ifTSl . 

Algorithm 1 Algorithm for minimum crossing contra-polymatroid intersection. 
1: Initially, set E' = E, F = (/},W = I, b[ = bi, for alH G / 
2: while E' ^(/}do 

3: Compute an optimal extreme point solution x* of the LP{E', F, W); 
4: for all e G withx*(e) = Odo 
5: E' ^E'\ {e} 
6: end for 

7: foralle G withx*(e) > ^do 
8: F ^ FU{e};E' ^ E'\{e} 
9: b'--^bi- x*{e), for allieW with e e Ei 
10: end for 

11: for aU i G with \Ei D E'\ < \2b'^ + A - 1 do 
12: W ^W\ {i} 

13: end for 
14: end while 

15: Return the incidence vector of F; 



Note that this algorithm rounds variables of value x*(e) > ^ to 1, and hence we loose a factor of two 
in the cost and in the degree bounds. Theorem [5] follows as a consequence if we can show that in each 
iteration, either some variable can be rounded, or some constraint can be dropped. 

Lemma 14 If x* G is an optimal extreme point solution to the above LP for crossing contra- 
polymatroid intersection, with < x*{e) < ^ for all e ^ E, then there exists i € W such that 

lEiHE'l < [26^1 + A - 1 

Proof: Let 7^' = {x{E' n S)\x*{S n E') = ri{S) - \S n F\, S ^ E} for i = 1,2 denote the tight 
sets from the first two constraints of the LP. Let B' = {x{E' n Ei)\x*{Ei n E') = b'^, i £ W} denote 
the tight degree constraints. Since x* is an extreme point solution (and < x* < 1), there exist linearly 
independent tight sets Ti C 7^', 7i C and B ^ B' such that \E'\ = \Ti\ + \T2\ + \B\. 

Since x* is modular and ri{S) — \S H F\ (for i = 1, 2) are supermodular on 2^, it can be assumed 
(again, using uncrossing arguments) that each of (71, ^) and (Ti, Q) forms a chaii^l. The following 
claim goes back to a similiai^ result for spanning trees as stated in IH. 

Claim 15 For each i = 1,2, we have \Ti\ < x*{E'); additionally if\Ti\ = x*{E') then E' G %■ 

Proof: We prove the claim for i = 1. Let 71 = {Si C ... C Sk} where Sk ^ E' . Let 5*0 = and 
consider an arbitrary pair of subsequent chain elements Si C S'j+i, for any i G {0, 1, . . . , A; — 1}. Since 

"A family (£, C) is a chain iff for every X,Y a C, either X C F or y C X. 
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xl>0 for all ee E'it follows that x*{Si+i \ Si) > 0. Hence, by the integrahty of ri{S) - |5 n F| and 
tight constraints Si and 5'j+i, 

x*{Si+i \ Si) = x*{Si+i) - x*{S^) = n{S^+i) - \Si+i nF\- rl{S^) + \s^ n F\ > 1. 

Summing over i = 0, . . . , A; — 1 we therefore obtain the inequality: 

x*{E') > x*{Sk) = \Si)>k = \Ti\, 

with equality only if E' = Sk- ■ 

We now proceed with the proof of Lemma [I4l Suppose (for a contradiction) that for all i £ W, 

\Ei n E'\ > \2b'^ + A. For each i G W, define Sp^ := EeGi?'n£;,(l " = \E' n Ei\ - 2x*{Ei). 
Then we have Sp^ > \E' n Ei\ - 2h\ > \E' n Ei\ - \2b'^ > A. Hence J2i(^w SPi > A • \W\. 

For each e G i?', let rg := |{i S : e G < A the maximum element frequency. Note also that 
< 1 - 2x* < 1 for each e G E'. Now, 

J^Sp, = ^re-(l-2x:)<A. ^(l-2x:) 

iGiy e£E' e&E' 

= A • (l^'l - 2 • x*{E')) < A • (l^'l - ITil - ITal) 

The last inequality uses Claim [15] Note that equality holds above only if £" G 71 n Ti (by Claim [T5]). 
which would contradict the linear independence of 71 and Ti- Thus we have: 

^ Spi < A • {\E'\ - \Ti\ - \T2\) = A • 1^1 < A • \W\. 

However this contradicts the assumption \E' r\Ei\ > [26^] + A for all i G W. ■ 

Proof: [TheoremlU Lemma[T4limplies that an improvement is possible in each iteration of Algorithm[T] 
Since we only round elements that the LP sets to value at least half, the cost guarantee is immediate. 
Consider any degree bound i G /; let 6- denote its residual bound when it is dropped, and F' (resp. E') 
the set of chosen (resp. unfixed) elements at that iteration. Again, rounding elements of fractional value 
at least half implies | £'j n F' | < [2bi — 2b'^\ = 2bi — [26^] . Furthermore, the number of Fj-elements in the 
support of the basic solution at the iteration (ie. E') when constraint i is dropped is at most [26^] + A — 1. 
Thus the number of Fj-elements chosen in the final solution is at most \Ei n F'\ + \Ei n E'\ < 26j — 
[26^1 + [26'.] + A- l = 2-6i + A-l. H 

Tight Example. We note that the natural iterative relaxation steps (used above) are insufficient to obtain 
a better approximation guarantee. Consider the special case of the crossing bipartite edge cover problem. 
The instance consists of graph G which is a 4n-length cycle, with its edges partitioned into two perfect 
matchings Ei and E2. There is a degree-bound of n on each of Ei and E2; so A = 1. Consider the 
fractional solution to the LP-relaxation that assigns value of ^ to all edges. It is indeed a fractional 
edge-cover since each vertex is covered to extent one. The degree-bounds are cleai^ly satisfied. It is also 
an extreme point: note that this is the unique fractional solution minimizing the all-ones cost vector. 
For this extreme point solution, the largest edge- value is i, and the support-size (i.e. 2n) of its degree- 
constraints is twice their bound (i.e. n). Thus the iterative relaxation must either pick a half-edge or drop 
a degree-constraint that is potentially violated by factor two. 
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5 Minimum Crossing Lattice Polyhedra 



Before formally defining the lattice polyhedra problem, we need to introduce some terminology. We 
use notation similar to llT4l . Let [F, <) be a partially ordered set with F ^ %. We consider a lattice 
{J^, <)^ where there are two commutative binary operations, meet A and join V, that are defined on all 
pairs A,BeF, such that: 

AaB < A,B < AVB 

Note that our definition is more general than the usual definition of a lattice, since the join Av B is 
not required to be the least common upper bound of A and B. A function r : F ^ Z+ is said to be 
supermodular on [F, <, A, V) iff: 

r{A) + r{B) < r{A A B) + r{A V B), for all A,B eF 

Given a supermodular function r : F ^ Z+, a ground set E, a cost function c : E ^ M+, and a 
set-valued function p : F ^ 2^ satisfying: 

1. Consecutive property: If ^ < S < C then p{A) n p{C) C p{B), 

2. Submodularity: For all A,B €F, p{A V B)LI p{A A B) Q p{A) U p{B), 
the lattice polyhedron problem is defined as the following integer program: 

minic'^-xl ^ > r(5), V5 G -F; x G {0, 1}^ 
I egp(5) 

Definition 16 (Minimum crossing lattice polyhedron) Given a lattice polyhedron {E, {F, <), r, p, c) 

as above, and lower/upper bounds {aj}ig/ and {6j}ie/ on a collection {Ei C E}i^j, the goal is to 
minimize: 

-x] ^ Xe> r{S), yS eF; < x{Ei) < hi, yi e I; xe {0, 1}^ 

eep(S) 

We already mentioned in the introduction that several discrete optimization problems fit into the 
lattice polyhedron model (see e.g. |[29l ). 

For example, in the contra-polymatroid intersection problem with two supermodular rank functions 
ri,r2 : 2^ — )• R, the lattice {F, < ) consists of two copies 5' and S" for each subset S Q E, with partial 
order: 

5' < T" and {S QT =^ S' < T', S" > T"); V 5, T C 2^. 

This is easily seen to satisfy the consecutivity and submodularity properties. The rank function r for the 
lattice polyhedron has r(S") = ri{S) and r{S") = r2{S), for all S C E. 

In the planar min-cut problem, recall that F consists of all s — t paths in the given s, t-planar graph 
G. The partial order sets for any pair oi s — t paths P, Q, 

P < Q <^=^ P "below" Q in the planai" representation. 

The induced lattice turns out to be consecutive and submodular. The rank function is the all-ones func- 
tion. For more details on the relation between planar min cut and lattice polyhedra, the reader is referred 

tons. 
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5.1 Integrality gap for general crossing lattice polyhedra 

We first show that there is a bad integrality gap for crossing lattice polyhedra. Consider the planar 
min-cut instance on graph G = {V,E) in Figure |4] with vertices s,t E y as shown. Define edge-sets 
Ei := {{vi-i,Uij)}j^j^ U{(^«> '"i,j)}^=i foi" s^ch z G [k]; here we set vq = s and Vk = t. There are only 
degree upper-bounds in this instance, namely bound of one on each {Ei}f^^. Note also that A = 1 in 
this instance, and size of the ground-set n = \E\ = 0(/c^). 




Figure 4: The integrality gap instance for crossing planar min-cut. 

Consider the LP solution that sets Xe = for every edge e & E. It is clearly feasible for the rank 
constraints (every s — t path has x-value one). Furthermore, x{Ei) = \Ei\/{2k) = 1 for all i ^ [k]; i.e. 
the degree constraints are also satisfied. Hence the LP relaxation is feasible. 

On the other hand, consider any integral solution I E that has \I r\ Ei\ < A; — 1 for alH G [k\. 
It can be checked directly that there is an s — t path using only edges E\ I. Thus any integral feasible 
solution J must have maxjgj;.] | J n E'j| > k, i.e. it violates some degree-bound by at least an additive 
k — I = Q.{^/n) term. 

5.2 Algorithm for crossing lattice polyhedra satisfying monotonicity 

Given this bad integrality gap for general crossing lattice polyhedra, we are interested special cases 
that admit good additive approximations. In this section we consider lattice polyhedra that satisfy the 
following monotonicity property, and provide an additive approximation. 

(*) S<T =^ \p{S)\ < \p{T)\, for all S,TeT 

As noted earlier, this property is satisfied by all matroids, and so our results generalize that of Kiraly 
et al. |[T9l . In the rest of this section we prove Theorem |6] The algorithm is again based on iterative 
relaxation. At each iteration, we maintain the following: 

• F C E of elements that have been chosen into the solution. 

• E' E \ F of undecided elements. 

• W CI [m] of degree bounds. 

Initially E' = E, F = $ and W = [m]. In a generic iteration with E' , F, W, we solve the following 
LP relaxation on variables {xe \ e G E'}, called CP\^[{E\ F, W): 

T 

mm c X 

xip{S)) > r{S) - \F n p{S)\, V5G-F 
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ai-\FnEi\ <x{Ei) <bi-\FnEi\, e W 

< Xe < 1, Ve G E'. 

Consider an optimal basic feasible solution x to the above LP relaxation. The algorithm does one of 
the following in iteration (£", F, W), until E' = W = ^. 

1. If there is e e E' with Xe = 0, then E' ^ E' \ {e}. 

2. If there is e G with Xe = 1, then F ^ F U {e} and E' ^ E' \ {e}. 

3. If there isi^W with \Ei n E'\ < 2A, then W ^W\ {i}. 

We note that this algorithm is a natural extension of the one for matroids iTTOl and the one for spanning 
trees lOTI . However the correctness proof (next subsection) relies only on properties of lattice polyhedra 
and the monotonicity property (*). 

5.3 Proof of Theorem g] 

Assuming that one of the steps ([T])-© applies at each iteration, it is clear that we obtain a final solution F* 
that has cost at most the optimal value, satisfies the rank constraints, and violates each degree constraint 
by at most an additive 2A — 1. We next show that one of ([T])-© applies at each iteration (£", F, W). 

Lemma 17 Suppose (F, < ) is a lattice satisfying the consecutive and submodular properties, and con- 
dition {*), function r is supennodular, and x is a basic feasible solution to CVi^f with < Xg < Ifor 
all e £ E' . Then there exists some i G W with {EiCi E'\ < 2 A. 

We first establish some standai^d uncrossing claims (Claim [18] and Lemma [19]). before proving this 
lemma. We also need some more definitions. Two elements A, B G are said to be comparable if either 
A < B or B < A; they are non-comparable otherwise. A subset £ C J^" is called a chain if C contains 
no pair of non-comparable elements. Note that a chain in F does not necessarily correspond to a chain 
in 2^ (with the usual subset relation) under mapping p. 

Let r'{S) := r{S) — \F H p{S)\ for all 5 G -F denote the right hand side of the rank constraints in 
the LP solved in a generic iteration (£", F, W). 

Claim 18 r' is supermodular. 

Proof: This follows from the consecutive and submodulai^ properties of lattice (F, <). Consider any 
A,B e F, and 

\Fr\pA\ + \Fr\pB\ = \Fr\{pAyj pB)\ + \Fr\{pAr\pB)\ 

> \Fn{pAAB^pAvB)\ + \Fn{pAnpB)\ 

> \Fn {pAAB u pavb)\ + \Fn {pAAB n pa\/b)\ 
= \FnpAAB\ + \FnpAvB\ 

The second inequality follows from submodularity (i.e. pa^Pb 5 PAab^PAvb), and the third inequality 
uses the consecutive property paab H pAyB ^ PA, Pb (since Af\B < A,B < A\/ B). This combined 
with supermodularity of r imphes r'{A) + r'{B) < r'{A A B) + r'{A V B) for all A,B e F. ■ 

For any element Ag F, let G {0, 1}-^' be the incidence vector of p{A) C E'. Let T := {A G F \ 
x{pa) = t'{A)} denote the elements in F that correspond to tight rank constraints in the LP solution x of 
this iteration. Using the fact that r' is supermodular (from above), and by standai^d uncrossing arguments, 
we obtain the following. 
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Lemma 19 IfS,T e F satisfy x{ps) = r'{S) and x{pt) = r'{T), then: 



x{p{S AT))=r'{S AT) and x{piS V T)) = r'{S V T) 

Moreover, xiS) + x(T) = xiS A T) + xiS V T). 
Proof: We have the following sequence of inequalities: 

r'{S AT) + r'{SVT) < x{psat) + x{psvt) 

= x{psAT n PS\jt) + x{psAT U Psvt) 

< x{psAT n PS\/t) + x{ps u pt) 

< x{ps n pt) + x{ps u pt) 
= x{ps) + x{pt) 

= r'{S)+r'{T) 

< r'{S AT) + r'{SVT) 

The first inequality is by feasibility of x, the third inequality is the submodular lattice property, the 
fourth inequahty is by consecutive property, and the last inequality is supermodularity of r'. Thus we 
have equality throughout, in particular x{p{S V T)) = r'{S V T) and x{p{S A T)) = r'{S A T). Finally 
since Xe > for all e G E', we also have x{S) + x{T) = x{S A T) + x{S V T). m 

Given Claim[T8]and Lemma[T9l we immediately obtain the following (see eg. |[29l . Chapter 60). 

Lemma 20 (129J) There exists a chain C T such that the vectors {x{A) I ^ £ ^} ^''^ linearly 
independent and span {xiB) \ B G T}. 

We are now ready for the proof of Lemma [T71 

Proof: [Lemma [13 \E'\ is the number of non-zero variables in basic feasible x. Hence there exist 
tight linearly independent constraints: C F corresponding to rank-constraints and B W degree- 
constraints, such that \E'\ = \C\ + \B\. Furthermore, by Lemma l20l£ is a chain in F, say consisting of 
the elements Si < S2 < • • • < Sk- We claim that, 

\p{Sj) \ (iJlZlpiSt)) I > 2, for each l<j<k (2) 

The above condition is clearly true for j = 1: since x{p{Si)) = r'{Si) > 1 (it is positive and integer- 
valued), and Xe < I for all e G E'. Consider any j > 2. By the consecutive property on St < •S'j-i < Sj 

(for any 1 < t < j-l), we have p{S,)np{St) C p{S,^i). So, p{Sj)\(u{zlp{St)) = p{Sj)\p{Sj^i). 
We now claim that \p{Sj) \ p{Sj^i)\ > 2, which would prove Since < Sj, assumption (*) 
implies that there is at least one element e G p{Sj) \ p{Sj^i). Moreover, if this is the only element, i.e., 
if p{Sj) \ p{Sj^i) = {e}, then p(5j_i) = p{Sj) \ {e} must be true (again by property (*)). But this 
causes a contradiction to the non-integrality of x^: 

xe = x {p{Sj)) - X (piSj^i)) = r' {p{S,)) - r' {p{S,^i)) G Z. 

Now, equation (|2j) implies that k = \C\ < ^-i-. Hence \E \ < 2\B\. 

Suppose (for contradiction) that \Ei n £"| > 2A + 1 for all i G W. Then Y^ieW ^ 
(2A + 1) • \W\. Since each element in E' appears in at most A sets {Ei}i^\y, we have A • \E'\ > 
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Y^idw \Ei ^ ^ (2^ + 1) • 1^1- Thus \E'\ > 2\W\ > 2\B\, which contradicts \E'\ < 2\B\ from 



We are now able to prove the main result of this section: 
Proof: [Theorem|6l Since the algorithm only picks 1-elements into the solution F, the guarantee on cost 
can be easily seen. As argued in Lemma [iTJ at each iteration (£", F, W) one of the Steps ([B-® apply. 
This implies that the quantity l-E'l + \ W\ decreases by 1 in each iteration; hence the algorithm terminates 
after at most \E\ + |/| iterations. To see the guarantee on degree violation, consider any z G / and let 
{E', F, W) denote the iteration in which it is dropped, i.e. Step ^ applies here with \Ei n E'\ < 2 A 
(note that there must be such an iteration, since finally W = 0). Since a degree bound is dropped at this 
iteration, we have < < 1 for all e G £" (otherwise one of the earlier steps ([Til or Q applies). 

1. Lower Bound: a, - |F n Ei\ < x{Ei n E') < n < 2A, i.e. Oj < |F n Ei\ + 2A - 1. The 
final solution contains at least all elements in F, so the degree lower bound on Ei is violated by at 
most 2A - 1. 

2. Upper Bound: The final solution contains at most \F r\ Ei\ + \E' r\ Ei\ elements from Ei. If 
EiDE' = 0, the upper bound on Ei is not violated. Else, < x{Ei D E') < bi - \F D Ei\, i.e. 
k > l + \FnEi\, and |F n Ei\ + \E' n Ei\ < bi + 2A - 1. So in either case, the final solution 
violates the upper bound on Ei by at most 2A — 1. 

Observing that all the steps ([Hi-© preserve the feasibility of the ^'Pjat' follows that the final solution 
satisfies all rank constraints (since E' = % finally). H 

5.4 Algorithm for inclusion-wise ordered lattice polyhedra 

We now consider a special case of minimum crossing lattice polyhedra where the lattice F is ordered 
by inclusion. I.e. the partial order in the lattice is the usual subset relation on 2^. This class of lattice 
polyhedra clearly satisfies the monotonicity property (*), so Theorem [6] applies. However in this case, 
we prove the following stronger guai^antee for the setting with only upper bounds. This improvement 
comes from the use of fractional tokens in the counting argument, as in |@1 (for spanning trees) and 1191 
(for matroids). 

Theorem 21 If the underlying lattice of the minimum crossing lattice polyhedron problem is ordered by 
inclusion and only upper bounds are given, then there is an algorithm that computes a solution of cost at 
most the optimal, where all rank constraints are satisfied, and each degree bound is violated by at most 
an additive A — 1. 

The algorithm remains the same as the one above for Theorem [6] In order to prove Theorem |2T] it 
suffices to show the following strengthening of Lemma [TTl 

Lemma 22 Suppose {F, < ) is a lattice satisfying condition 



above. 



S<T 




function r is supermodular, and x is a basic feasible solution to ^Fi^f with < Xg < ^for all e G E'. 
Then there exists some i £ W with \Ein E'\ < 6^ + A — 1. 
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Proof: The proof is very similiai" to the proof of Lemma [141 Cleai^ly, since T is ordered by inclusion, 
the consecutivity and submodularity property are satisfied. Since x is a basic feasible solution, there exist 
linearly independent tight rank function- and degree bound constraints T and B '^W such that 

\E'\ = \T\ + \B\. 

Using uncrossing arguments, we can assume that (T, < ) forms a chain 

T = {Ti<T2<...< Tk}. 

Consider an arbitrary pair Ti < Tj+i in T, where i G {1, . . . ,k — 1}. Since Xg > for all e £ E and 
p{Ti) C p{Ti+i), it follows that < x{p{Ti+i) \ p{Ti)) and therefore, by the integrality of r, 

x(p(r,+i) \ p(r,)) = x{p{T,+^)) - x{p{T,)) = r{T,+,) - r(r,) > 1. 

By a similar argument, x{p{Ti)) > 1. Thus, 

k-l 

x{E) > x{p{n)) = ^(pm+i) \ p{Ti)) + x(p(ri)) >k = \T\ 
1=1 

with equality only if E = p{Tk). This implies that 

\E'\-x{E) = \T\ + \B\-x{E) <\B\. (3) 
Let E[ = E' n Ei. To prove the statement of the Lemma, it suffices to show: 

Y,{\E[\-h[) = Y,{\E[\-x{E,))<m\- 

In order to prove this, define Ae = |{i S TV | e G Ei]\ and consider the derivations 

Y,m-x{E,)) = ^ 5](l-x,) = 5]Ae(l-x,) 

= Aj^a-Xe)- J](A-Ae)(l-Xe) 
^ A|e| - J^(A- Ae)(l-Xe) 

e^ 

= A\W\ - A\W \B\- ^(A - Ae)(l - xe) < A\W\. 

eeE 

Note that equality can only hold if E = p{Tk) and A|Ty \ B\ + Y^eaEi^ " ^e)(l - x^) = 0. The latter 
can only be true if = \ W\ and Ag = A for each e £ E. But this would imply that 

5]x^' = Ax^ = Ax^^ 

where S {0, 1}-^^^ is the incidence vector of S £ T with xf = 1 iff e S p{S). However, this 
contradicts the fact that the constraints T and B are linearly independent. ■ 



Acknowledgement: We thank Mohit Singh 11301 for the integrality gap for general MCST, and Chandra 
Chekuri for finding an error in the arborescence result in an earlier version Q of this paper. 



22 



References 

[1] B. Anthony, V. Goyal, A. Gupta, and V. Nagarajan, A Plant Location Guide for the Unsure: Ap- 
proximation Algorithms for Min-Max Location Problems, Math. Oper. Res., 35, 2010, 79-101. 

[2] S. Arora, L. Babai, J. Stem, and Z. Sweedyk, The hardness of approximate optima in lattices, codes, 
and systems of linear equations , J. Comput. Syst. Sci., 54(2), 1997, 317-331. 

[3] V. Arya, N. Garg, R. Khandekar, A. Meyerson, and K. Munagala, and V. Pandit, Local search 
heuristics for k-median and facility location problems, SIAM J. on Computing, 33(3), 544-562, 
2004. 

[4] N.Bansal, R. Khandekar and V. Nagarajan, Additive Guarantees for Degree-Bounded Directed Net- 
work Design, SIAM J. Comput., 39(4), 2009, 1413-1431. 

[5] N.Bansal, R. Khandekar, J. Konemann, V. Nagarajan, and B. Peis, On Generalizations of Network 
Design Problems with Degree Bounds, In IPCO, 2010. 

[6] V. Bilo, V. Goyal, R. Ravi and, M. Singh, On the crossing spanning tree problem. In APPROX 
2004, 51-60. 

[7] M. Charikar, S. Guha, E. Tardos, and D.B. Shmoys, A Constant-Factor Approximation Algorithm 
for the k-Median Problem, J. Comp. and Syst. Sciences, 65(1), 129-149, 2002. 

[8] K. Chaudhuri, S. Rao, S. Riesenfeld, and K. Talwar, What would Edmonds do? Augmenting paths 
and witnesses for degree-bounded MSTs, Algorithmica, 55(1), 2009, 157-189. 

[9] K. Chaudhuri, S. Rao, S. Riesenfeld, and K. Talwar, A push-relabel approximation algorithm for 
approximating the minimum-degree MST problem and its generalization to matroids, Theor. Com- 
put. Sci., 2009, 410(44), 4489-4503. 

[10] B. Chazelle, The Discrepancy Method: Randomness and Complexity, Cambridge University Press, 
2000. 

[11] C. Chekuri, J. Vondrak and Rico Zenklusen, Dependent Randomized Rounding for Matroid Poly- 
topes and Applications, [http://arxiv.org/abs/090974348l 2009. 

[12] U.Faigle and B. Peis, Two-phase greedy algorithms for some classes of combinatorial linear pro- 
grams. In SODA, 2008, 161-166. 

[13] U. Faigle, W. Kern and B. Peis, A ranking model for cooperative games, convexity and the greedy 
algorithm, working paper. 

[14] A. Frank, Increasing the rooted connectivity of a digraph by one. Math. Programming 84 (1999), 
565-576. 

[15] M.X. Goemans, Minimum Bounded-Degree Spanning Trees, In FOCS, 273-282, 2006. 

[16] F. Grandoni, R. Ravi, M. Singh, Iterative Rounding for Multiobjective Optimization Problems, In 
ESA, 2009, 95-106. 

[17] A. Hoffman and D.E. Schwartz, On lattice polyhedra. In Proceedings of Fifth Hungarian Combi- 
natorial Coll. (A. Hajnal and V.T. Sos, eds.), North-Holland, Amsterdam, 1978, pp. 593-598. 



23 



[18] K. Jain, A factor 2 approximation algorithm for the generalized Steiner network problem, Combi- 
natorica, 2001, 39-61. 

[19] T. Kiraly, L.C. Lau and M. Singh, Degree bounded matroids and submodular flows. In IPCO 2008, 
259-272. 

[20] P. N. Klein, R. Krishnan, B. Raghavachaii and R. Ravi, Approximation algorithms for finding low 
degree subgraphs, Networks, 44(3), 2004, 203-215. 

[21] J. Konemann and R. Ravi, A matter of degree: Improved approximation algorithms for degree 
bounded minimum spanning trees, SIAM J. on Computing, 31:1783-1793, 2002. 

[22] J. Konemann and R. Ravi, Primal-Dual meets local search: approximating MSTs with nonuniform 
degree bounds, SIAM J. on Computing, 34(3):763-773, 2005. 

[23] B. Korte and J. Vygen. Combinatorial Optimization. Springer, New York, 4th ed., 2008. 

[24] L.C. Lau, J. Naor, M. R. Salavatipour and M. Singh, Survivable network design with degree or 
order constraints, SIAM J. on Computing, 39(3), 2009, 1062-1087. 

[25] L.C. Lau and M. Singh, Additive Approximation for Bounded Degree Survivable Network Design, 
In STOC, 2008, 759-768. 

[26] Z. Nutov, Approximating Directed Weighted-Degree Constrained Networks. In APPROX, 2008, 
219-232. 

[27] R. Ravi, M.V. Marathe, S.S. Ravi, D.J. Rosenkrantz, and H.B. Hunt, Many birds with one stone: 
Multi-objective approximation algorithms. In STOC, 1993, 438-447. 

[28] R. Ravi and M. Singh, Delegate and Conquer: An LP-based approximation algorithm for Minimum 
Degree MSTs. In ICALR 2006, 169-180. 

[29] A. Schrijver, Combinatorial Optimization, Springer, 2003. 

[30] M. Singh, Personal Communication, 2008. 

[31] M. Singh and L.C. Lau, Approximating minimum bounded degree spanning trees to within one of 
optimal, In STOC, 2007, 661-670. 



24 



